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We investigate superresonant scattering of acoustic disturbances from a rotating acoustic black 
hole in the low frequency range. We derive an expression for the reflection coefficient, exhibiting its 
frequency dependence in this regime. 



I. INTRODUCTION 



' Progress in understanding black hole radiance, both in the semiclassical and quantum gravitational regimes, have 

. suffered greatly in the past, due to the lack of experimental feedback. Black hole radiance consists of spontaneous 
radiation (Hawking radiation) [l]as well as stimulated emission (superradiance) [2], [3], [4], [5]. For astrophysical black 
holes, the Hawking temperature is invariably smaller than the Cosmic Microwave Background temperature. So the 
black hole mostly accretes rather than radiates spontaneously. Similarly, superradiant emission is also hard to observe 
from within a background of x-rays being emitted from accretion processes. Due to these reasons physicists have 
been looking for alternative physical models of gravity so that experiments can be devised for indirect verification 
of phenomena that are predicted by the semiclassical theory of gravity. These alternative physical models are called 

. analog models of gravity. 

Among the first of such analog gravity models, Unruh's [6] was the pioneering idea of analog models based on fluid 
mechanics. The idea of such models arose from a remarkable observation regarding sound wave propagation in a fluid 

' in motion. Unruh showed that if the fluid is barotropic and inviscid, and the flow is irrotational, the equation of 
t~H , motion of the acoustic fluctuations of the velocity potential is identical to that of a minimally coupled massless scalar 

■ field in an effective (3+1) dimensional curved spacetime, and is given by, 

m ; 

£o : = -L a M (v=5 <T d v ) * = o (i) 
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where g^ u is the 'acoustic' metric and is the fluctuation of the velocity potential. The most remarkable fact 
' about this effective spacetime 'seen' by the acoustic fluctuations is that the acoustic metric is curved and Lorentzian in 
signature, even though the fluid flow is governed by non-relativistic physics in a flat Minkowski background spacetime. 
This particular feature of these models raises the possibility of being able to mimic a gravitational system in its 
kinematical sector. Since the fluid has barotropic equation of state, the pressure does not generate vorticity. So an 
• *h , initially irrotational flow will remain irrotational. Viscosity of the fluid has been ignored in order to maintain the 
/\i ' Lorentzian character of the effective acoustic geometry. 

In acoustic analog models of gravity, the acoustic geometry entirely depends on the the background fluid motion. 
So the changes of the acoustic geometry are also related to that of the background fluid motion. Since the background 
fluid motion is governed by the Euler equation and the continuity equation, not by Einstein's equation, it may mimic 
only the kinematical aspects of gravity, and not its dynamical ones. If the flow contains a well-defined region where 
it is supersonic, that region resembles an ergoregion of a rotating black hole, characterized by the time translation 
Killing vector turning spacelike. Within this region, there may exist a surface where the inward normal velocity of the 
fluid exceeds the local speed of sound at every point. Such a surface therefore traps all acoustic disturbances within 
it, resembling an outer trapped surface, and can therefore serve as the sonic or acoustic horizon of the analogue black 
hole. 

Unruh [6] has used this idea of an acoustic analogue black hole horizon, to show that quantized acoustic perturbations 
(phonons) may be emitted from it precisely like the emission of photons from a black hole horizon via Hawking 
radiation. The Hawking temperature is likewise given by the surface gravity at the sonic horizon. This is acceptable 
provided we continue to ignore the backreaction of the phonon radiation on the background metric, which is a 
legitimate thing to do in a fixed-metric problem. Unruh's work has led, in the last couple of years, to several 
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proposals for the experimental verification of acoustic Hawking radiation. The Hawking temperature for typical 
acoustic analogues has been estimated to be in the nano-Kelvin range; with improved experimental techniques and 
equipment, this is not outside the realm of realistic possibilities in the near future. 

One slightly unsavory feature of the Unruh approach has to do with the necessity of quantizing linearized acoustic 
disturbances and the manner in which they are quantized. In a classical fluid, there is no compelling physical reason 
to quantize acoustic disturbances. So quantizing acoustic perturbations just to demonstrate Hawking radiation from 
the horizon of acoustic black holes appears somewhat artificial. In a quantum fluid like superfluid helium on the other 
hand, there are phonon excitations of the fluid background itself. On top of this one may have acoustic perturbations 
which are quantized in terms of phonons. However, in that case, the division between the background and linear 
acoustic fluctuations tends to get a bit hazy. 

There is however an alternative phenomenon which does not necessitate quantization of acoustic disturbances for it 
to occur, although a full quantum description is possible in principle. This is the phenomenon called Superresonance 
by us [7] which occurs in a curved Lorentzian acoustic geometry if the acoustic spacetime is that of a rotating black 
hole. Recall that such a black hole spacetime admits a region - the ergoregion - containing the event horizon, where 
the time translation Killing vector turns spacelike. The existence of this region is crucial for the Penrose process of 
energy extraction from the black hole, at the cost of its rotational energy. Superresonance is simply an acoustic wave 
version of the Penrose process, wherein a plane wave solution of a massless scalar field in the black hole background 
is scattered from the ergoregion with an amplification of its amplitude. The energy gained by the wave is at the cost 
of the rotational energy of the black hole. 

In the earlier paper [7], superresonance has been shown to occur in a certain class of analogue 2+1 dimensional 
rotating black holes. Linear acoustic perturbations in such a background are shown to scatter from the ergoregion 
with an enhancement in amplitude, for a restricted range of frequencies of the incoming wave. However, while the 
plausibility of superresonance was established in that earlier work through the analysis of the Wronskian of the radial 
equation of the acoustic perturbations, the actual frequency dependence of the amplification factor was not worked 
out. In this paper this gap is filled. We present a more detailed quantitative analysis of this phenomenon, as also 
a somewhat different way of arriving at the results of the previous paper. An explicit expression is derived for the 
reflection coefficient as a function of the frequency of the incoming wave for a certain low frequency regime. This 
expression is expected to be useful for possible future experimental endeavors to observe superresonance. 

The plan of the paper is as follows: in section II, certain aspects of the acoustic black hole analogue we work with 
are reviewed. The demonstration of superresonance is repeated with a choice of coordinates different from those used 
in the earlier work [7]. In section III, detailed quantitative analysis of the reflection coefficient is presented, explicitly 
exhibiting its frequency dependence, albeit within a certain restricted low frequency regime. We conclude in section 
IV with a sketch of the outlook on this class of problems. 



II. DRAINING VORTEX FLOW AND SUPERRESONANCE 



In order to demonstrate superresonance, a certain choice is to be made of the velocity potential of the fluid, which 
functions as a suitable acoustic analogue of a Kerr black hole. In other words, the flow must have a well-defined 
(ergo)region of transonic flow containing a sonic horizon as discussed in the Introduction. In 2+1 dimensions, a 
velocity profile with these properties is the 'draining vortex' flow with a sink at the origin, given by [8] 



_ A . B 

v = r H 

r r 



(2) 



where A, B are real and positive and (r, <j>) are plane polar coordinates. The metric corresponding to this effective 
geometry is , 



ds 



A 2 + B 2 \ , 2 2 A 
1 dt 2 + 



dr dt - 2B d6 dt + dr 2 + r 2 d6 2 



(3) 



A two surface, at 



in this flow, on which the fluid velocity is everywhere inward pointing and the radial 



component of the fluid velocity exceeds the local sound velocity everywhere, behaves as an outer trapped surface 
in this acoustic geometry. This surface can be identified with the future event horizon of the black hole. Since the 
fluid velocity (2) is always inward pointing the linearized fluctuations originating in the region bounded by the sonic 
horizon cannot cross this boundary. As for the Kerr black hole in general relativ ity, the r adius of the boundary of 
the ergosphere of the acoustic black hole is given by vanishing of g oa , i.c, r e = \J A 2 + B 2 /c. If we assume that the 
background density of the fluid is constant, it automatically implies that background pressure and the local speed of 
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sound are also constant. Thus we can ignore the position independent pre- factor in the metric because it will not 
effect the equation of motion of fluctuations of the velocity potential. 

From the components of the draining vortex metric it is clear that the (2+l)-dimensional curved spacetime possesses 
isometries that correspond to time translations and rotations on the plane. The solution of the masslcss Klein Gordon 
equation can therefore be written as, 



$>(t,r,<p) = R(r) e 

where uu and m are real and positive. In order to make ^(t,r,(j)) single valued , m should take integer values. 
Then the radial function R(r) satisfies 

d 2 R(r) „ , . dR(r) „ , N „, , „ 
+ Pi(r) —r 1 + Qi(r) R(r) = 
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where, 



and 



dr 2 
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Observe that the equ.(5) is different from that used in [7] (equ. 8), since the diffeomorphisms given in equ. (5) of [7] 
have not been effected in the present case. 

Now we introduce tortoise coordinate r* through the equation, 



d 

dr* 



1 



A 2 



which implies that, 



A , 
r + — log 
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This tortoise coordinate spans the entire real line as opposed to r which spans only the half-line. The horizon at r = ^ 
maps to r* — > — oo , while r — > oo corresponds to r* — ► +oo . Let us now define a new radial function G(r) as, 



/ x H(r) 
R(r) = exp 
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Substituting this in equation (5) we observe that G{r) satisfies the differential equation 
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Now in terms of tortoise coordinate one obtains the modified differential equation as, 



d 2 H(r*) 



dr 



*2 



B m 



A 2 



in 



+ 



A 2 



H{r*) = 



(8) 

H(r) = 
(9) 

(10) 



We analyze this differential equation in two distinct radial regions, viz., near the sonic horizon, i.e., at r* — > — oo 
and at asymptopia, i.e., at r* — ► +oo . In the asymptotic region, the above differential equation can be written 
approximately as, 



+ 4 "<>-> = » 

dr* 2 c 2 



(11) 
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This can be solved trivially, 



H(r*) = R um exp (i - r*) + cxp (-i - r*) (12) 

The first term in the above equation corresponds to reflected wave and the second term to the incident wave, so that 
R is the reflection coefficient in the sense of potential scattering. Similarly , near the horizon the above differential 
equation can be written approximately as, 

qp+t^i^^ (13) 

where is the angular velocity of the sonic horizon. We impose the physical boundary condition that of the two 
solutions of this equation, only the ingoing one is physical, so that one has, 

tf(r*) = T_ex P [-i r *j (14 ) 

Now using these approximate solutions of the above differential equation together with their complex conjugates 
and recalling the fact two linearly independent solutions of this differential equation (10) must lead to a constant 
Wronskian, it is easy to show that, 

l-|^ m | 2 =(l-^) |T wro | 2 (15) 

Here R^ m and T wm are amplitudes of the reflection and transmission coefficient of the scattered wave respectively. 
It is obvious from the above equation that, for frequencies in the range < ui < m fin, the reflection coefficient 
has a magnitude larger than unity. This is precisely the amplification relation that emerges in our earlier analysis 
of superresonance [7] . The demonstration that the final physical result remains the same, despite using a different 
system of coordinates in the acoustic spacetime, implies that the analogue spacetime indeed exhibits the kind of general 
coordinate invariance associated with standard general relativity. This relation between reflection and transmission 
coefficients, however, only establishes the plausibility of superresonance for draining vortex flows. It is not sufficiently 
detailed to link up with possible forthcoming experimental observations. For that purpose, we need to know the 
frequency dependence of both the reflection and the transmission coefficients. This requires that we have to solve the 
differential equation explicitly. In the next section we analyze the radial differential equation in the low frequency 
limit u A/c 2 < 1. 



III. THE REFLECTION COEFFICIENT 

We begin with defining a new function L(r) through the relation 

m B 



R( r ) = -j exp 



i I u A , 

— ^g 



2 2 

r z c 
~Ai~ 
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We next introduce a new variable x = 1 . Then L(x) satisfies the differential equation 



d 2 L(x) 
dx 2 



x x + 1 



where, 



dL(x) 
dx 



X X + 1 C 4 



L(x) 
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(17) 
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In order to find a solution of this differential equation, we adapt a matching procedure employed by Starobinsky [3]; 
for the region 



4 



uj A 



x <C m 



the above differential equation can be approximately written as, 



d 2 L(x) | (1 
dx 2 



1 



x + 1 



dL{x) 
dx 
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91 

x 



1 



X + 1 



< 1 



(1-S 2 ) 



L(x) 



4x(x + 1) 



= 



(18) 



This is the standard Riemann-Papparitz equation with regular singular points at 0,-1, and at oo. We can convert 
this differential equation into the Hypergeometric form by the following substitution, 



L(x) = x a (x + if G(x) 
such that G(x) satisfies the Hypergeometric differential equation given by 

3?(1 - + I 7 - {a + + l)x ] - a0G(x) = 



dx 



where 



a 



. Q 



S .0 R S .Q 
2 2 2 2 
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(19) 
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The most general solution of this differential equation is , 

G{x)=C x 2F 1 {a,f3; Tl x)+C 2 x 1 ^ 2 F 1 (a + 1 - 7 , + 1 - T , 2 - T, x) 
where G\ and Ci are arbitrary constants to be determined. Near the horizon, x ~ 

G(x) — ► d x a ' + C 2 x~ a ' . 



(21) 



(22) 



Here we have absorbed (— l) 1 1 in C% ■ The first term represents ingoing wave and the second term represents 
outgoing wave near the horizon. Since nothing can come out of the horizon, the boundary condition at the horizon 
implies that C 2 = 0. Therefore the solution of equ.(18) reduces to 



L{x)=Cix a {x + lf 2 F 1 (a,P; Tl -x) 



(23) 



To use the matching procedure we need to know the behavior of this solution in the asymptotic region, i.c, 
at x — ► oo . This can be obtained by the transformation (x — > 1/x) for the Hypergeometric function in 
above equation. Using this rule we obtain, 



L{x) = Cix a (x + if 



r( 7 ) rpg - a) 

r(/3) r( 7 -a)' 



x~ a 2 Fi [a , a + 1 - 7 ; a + 1 - 



(-1) 



M-ff 2Fl (0,0 + l-,.,0 + l- a .,tH 



T(a) r( 7 -/3) 
So in the region x 3> Max( m 2 ,Q 2 ), we have 

L(x) — ► C lX a ' (x + 1) " 



_r(/3) r( 7 -a) 
Hence the approximate form of L(x) in the region, 

Max( m 2 ,Q 2 ) < x < m. 



r(7) r(/3- a ) x _ a + r( 7 ) r(a - 0) x _ 



r(a) r( 7 -/3) 



(24) 



(25) 
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is given by 



L(x) = 



Ci 



r( 7 ) T(f3 - a) s 

r(/3) r( 7 - a) 



r( 7 ) T{a-p) _s 
T{a) r( 7 -/3) 



where, 



Di + D 2 x 2 



Di = r( 7 )rqg-a) Ci 



(26) 



ft = 



r( 7 -a) 

r( 7 ) r(a - g) 
r(a) r( 7 -/3) 



The solution (26) is in a form where it is not obvious which combinations of the complex coefficients D\ , D 2 corre- 
spond to the ingoing and outgoing modes, This information can be retrieved from the energy flux per unit time radially 

across an arbitrarily chosen surface at a constant radial coordinate r in the region Max(m 2 1 Q 2 ) « i « (^c^~)- 

This is obtained by the surface integral 



FL 



= - ! * F k^^f—gd^ 
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(27) 



where is the energy momentum tensor for ^, is the normal to a constant r surface , fc M is the time translational 
killing vector and g is the determinant of the acoustic metric. 

Given that the outgoing mode of the wave corresponds to a positive value of radial component of the flux, and the 
ingoing mode just the opposite, one concludes that the amplitudes of the incident and reflected wave are proportional 
to [D\ — i D 2 ) and {D\ + i D 2 ) respectively. 

One can now write down an explicit expression for the amplification factor(AF); this is given by, 



AF = 1 — \R 



= 1 - 
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(28) 



ft follows from this that the reflection coefficient exceeds unity in the frequency range < u> < m £Ih ■ Comparing 
this with the equation (15), we can identify T wm in this frequency range as, 



IT r - 



2 w A 



\Ci\ 



s 



IDi -i D 2 



(29) 



Equations (28-29) exhibit the frequency dependence of both the reflection and transmission coefficients and hence 
of the amplification factor in the low frequency regime. If vortex motion with a sink is realized in the laboratory, our 
result can be used to provide an order of magnitude estimate for the amplification factor. According to our assumption 
these relations are valid in the frequency range where "C 1 , i.e, when wave length of the sound wave is much 
larger than the radius of the horizon. Clearly, this frequency range does not cover the entire superresonant range 
to to Q,h . It is interesting and of importance to find the behavior of the reflection coefficient near the critical point 
around u> = m Q,h where the reflection coefficient is expected to cross unity. We hope to report on this in the near 
future. 
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IV. OUTLOOK 



The link with possible experimental observation of superresonance is yet incomplete, despite the explicit expres- 
sions we have derived for the reflection and transmission coefficients. One needs to input into these results typical 
characteristics of draining vortex flows in fluids like supcrfluid helium which obey our assumptions of barotropicity, 
irrotationality and zero viscosity. These characteristics are to be derived from the macroscopic quantum mechanics 
of supcrfluid helium, and novel features like flux quantization are expected to play a non-trivial role. In particular, 
flux quantization may well lead to the phenomenon of discrete amplification discussed in [7], which would perhaps 
facilitate observation. We should mention that ergoregions in supcrfluid helium have been considered in ref. [9], [10], 
although from a somewhat different standpoint. These authors do not consider vortex flows with a drain, and as such, 
their acoustic spacetime does not share all characteristics of black holes, like an outer trapped surface. 

Another feature which we have ignored in the foregoing analysis is the existence of shock fronts at the interface 
of normal and transonic flows. Indeed, in astrophysical black holes, analysis of these shock fronts constitutes an 
important component of models used in observational work on these objects. The existence of these shock fronts 
are in part responsible for the non-observation of superradiant scattering in such black holes. The issue of shocks in 
acoustic black hole analogues thus assumes special significance for future work on this topic. 

One of us (PM) thanks J. Bhattacharjee, G. Bhattacharya and A. Chatterjee for helpful discussions. 
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